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ABSTRACT 

The purpose of this paper is to demonstrate concretely the role that distributed 
parameter models can play in CSI, in particular in combined structures-controls optimiza- 
tion problems of importance in preliminaiy design. Closed form solutions can be obtained 
for performance criteria such as rms attitude error, making possible analytical solutions of 
die optimization problem. This is in contrast to the need for numerical computer solution 
involving the inversion of large matrices in traditional finite element model use. Another 
advantage of the analytic solution is that it can provide much needed “insight” into 
phenomena that can otherwise be obscured — or difficult to discern from numerical 
computer results. 

As a compromise in level of complexity between a toy laboratory model and a real 
space structure we have chosen the lattice truss used in the EPS (Earth Pointing Satellite). 
The optimization problem chosen is a generic one: of minimizing the structure mass 
subject to a specified stability margin and to a specified upper bound on the rms attitude 
error (“tip response”), using co-located controller and sensors. Standard FEM treating each 
bar as a truss element is used, while the continuum model is anisotropic Timoshenko 
beam model. Performance criteria are derived for either model, except that for the 
distributed parameter model we obtain explicit closed form solutions. Numerical results 
obtained by the two models show complete agreement. Based on the continuum model 
we obtain a solution to the problem of optimal placement of actuators to minimize mean 
square attitude error. A canonical optimization problem is examined and shown to be 
trivial, and even capable of analytical solution, using the continuum model performance 
criteria formulas in contrast to the complex computer solutions based on FEM or 
truncated modal models currently in vogue. 
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Introduction 


The most voiced criticism against the use of continuum models for structures is that 
they are (a) impossible to derive for a realistic structure and (b) even if it could be done, 
calculations using the model are equally impossible. We shall show that both statements 
are false — at least in the CSI optimization problem — in particular in preliminary 
design. 

In combined controls-structures optimization, the optimization is the least difficult 
the real challenge is to derive expressions for the chosen performance criteria in terms of 
the controls/structures parameters. We shall show that such formulas are much simpler 
when continuum models are used — moreover in many cases we can derive explicit 
closed form expressions in terms of elementary function which can actually trivialize the 
optimization problem. In particular the techinques of optimization need no longer 
dominate. 

The purpose of this paper is to demonstrate concretely the role that distributed 
parameter models can play in CSI, in particular in combined structures-controls optimiza- 
tion problems of importance in preliminary design. Closed form solutions can be obtained 
for performance criteria such as rms attitude error, making possible analytical solutions of 
the optimization problem. This is in contrast to the need for numerical computer solution 
involving the inversion of large matrices in traditional finite element model use. Another 
advantage of the analytic solution is that it can provide much needed “insight" into 
phenomena that can otherwise be obscured — or difficult to discern from numerical 
computer results. 

As a compromise in level of complexity between a toy laboratory model and a real 
space structure we have chosen the lattice truss used in the EPS (Earth Pointing Satellite). 
This is described in Section 1. The optimization problem chosen is a generic one: of 
minimizing the structure mass subject to a specified stability margin and to a specified 
upper bound on the rms attitude error (“tip response”). The mathematical statement of 
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the performance criteria is given in Section 2. The first step is to evaluate the performance 
criterion for a given control configuration — we consider co-located sensor/controls only. 
The finite element model is described in Section 3, and the continuum model in Section 4. 
The dynamic state space model is seen to be identical in both cases except for state space 
dimension. Section 5 derives the performance criteria valid for either model, except that 
for the distributed parameter model we obtain explicit closed form solutions. Section 6 
compares the numerical results obtained by the two models, showing complete agreement. 
As a byproduct of our analysis, we obtain a solution to the problem of optimal placement 
of actuators to minimize mean square attitude error — in Section 7. An optimization 
problem per se — a canonical one — is treated in Section 8 and by virtue of our explicit 
formulas for performance indices in terms of structure/control parameter, shown to be 
trivial and even capable of analytical solution — in contrast to computer solutions using 
FEM or truncated modal models as in [9, 10]. 

We should note that structural engineers (Noor, et al [1, 2] and Renton [3]) have 
already voiced the advantages of continuum models in preliminary structure design — 

what is new here is the application to control design, to the Controls-Structures Interaction 
problem. 
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1. The Physical Article 

The physical structure (Figure 1) is a lattice of rectangular bays, each single-laced 
single-bay. Offset at each end is an antenna. The controllers are force and moment 
actuators with co-located attitude as well as rate sensors stationed at arbitrary locations 
along the structure. Table 1 is a breakdown of the parameters describing each bay. 


TABLE 1 


Element Properties 



Battens 

Longitudinal 

Bars 

Diagonal 

Bars 

Cross Bracing 
in Battens 

Length L 

b 

i 

d 

5 

Sectional 
Area A 



A d 

^8 

Elastic 
Modulus E 

E b 

El 

E d 

E S 

Mass Density m 

P b 

Pt 

Pd 

Pe 

Element 
Mass = p AL 

m b 


m d 


Element 

Stiffness = EAIL 

s b 

s t 

s d 

5 6 


The beam geometry is shown in Figure 2. By the “nominal structure, we shall mean the 
following choice of structural parameters: 

b = « ; d = 6 = V2 £ 

A b = A t = A d = Ag = A 

= E b = E t = E d = = E 

- Pi - Pt = Pd = P6 “ P 
L - nZ ; 


Zb 7 



Nominal values: 


n = 9 
£ = 3 m 
p = 3250 
E = 2.759 x 10 n 
/I - 2.468 X 10"V . 
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2. Performance Criteria 


As in previous CIS optimization studies (see [9] and the references therein) the 
objective is to minimize the total mass of the structure, including the controller mass, 
subject to meeting specified performance requirements, which we shall now describe in 

mathematical terms. 

The performance criteria chosen are: 

(i) the mean-square attitude error due to sensor noise (using co-located sensors/ 
actuators). As we shall see, this actually depends on the steady state “tip” response to step 
inputs so that “noise” notions can be eschewed if necessary. 

(ii) the “stability margin”: defined as the sum of the absolute values of the real 
parts of the closed-loop eigenvalues. This is one measure of stability among veiy many 
(see, e.g., [10]). We choose this one because it is essentially equivalent to any other one 
but has the advantage that we can derive a simple closed-form expression for it. 

Let us now define the criteria more precisely. First of all we assume the control 
law to be PD (“proportional plus derivative”) as in classical servo design. Let v p (f) denote 
the “displacement” or “attitude” vectors at the sensor locations and let v r (f) the rate. 

Then U{t) the control is defined to be 

£7(0 = av p (0 + Yv r ( 0 (2 ' 

where a and y are (positive) scalar “gains.” This is not altogether a “simplifying” 
assumption — that the scalar rate feedback is actually optimal is shown in [4], The 
beam-axis being the x } -axis, we have, with L denoting the beam-length. 

0 < x\ < L . 

Let /(0),/(L) denote the 6-DOF displacement vectors at the ends. Let 


U{ 0) 


u{L) 

v(0) 


v(L) 

w(0) 

; m = 

w(L) 

<M°) 


<M L ) 

<t>2 (0) 


<t>2 (Z-) 

<S> 3 (0) 
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where u, v, w are the 3-DOF displacements, 4>, the torsion angle about the beam axis, 

and 4 > 2 , 4> 3 about the mutually perpendicular axes. Then the mean square attitude error is 
defined by 

~ w(0) z + u(L ) 2 + v(0) 2 + v(L) 2 + w(0) 2 + w(L) 2 

, ^ ( 23 ) 

+ M 4>i (0) 2 + \r L \ 2 <t>, (L) z , 

where the bars denote time averages, |r| being the length of the moment arm as required. 

Under our feedback law (often referred to as “positive-definite” feedback) the closed- 
loop system is guaranteed to be stable. Assuming no damping in th structure (as we shall 
indeed do), the real parts of the closed-loop eigenvalues are guaranteed negative (see 
Section 5) — or if we assume the structure is already damped we have stability enhance 
ment corresponding to the increment, the real parts being now more negative. Let o, 
denote the real part increment corresponding to the rth closed-loop mode. Then the 
infinite series 

oo 

£ -o. (2.4) 

converges. Denote the sum by o,. We shall take this to be the “stability index” — the 
higher the index, the more stable. 
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3. The Finite Element Model 

Since most of the techniques in developing the FEM are standard, we shall only 
present the relevant numerical data. Each bar is taken as a truss element with 6 DOF. 
There are (13 + 5) elements per bay, and hence (13x9 + 5) = 122 elements for 9 bays. 
There are 40 nodes with 3 DOF each, so that the stiffness matrix A and the mass/inertia 
matrix M are 120x120. The state vector is thus 120x1. The displacements along the axis 


of the truss are then expressed (or, rather, extrapolated): 

u(kZ.0-b/2) + u(kZ,0+bl2) + u(kZ,-b!2,0 ) + u(kZ, +b/2, Q1 
u(kZ) = u(kZ, 0,0) = - L — 1 L 4 

\(kZ.Q-b!2) + v(kZ,0+b/2) + vfH. -fr/2.0) + v(*£, +6/2,0) 
v(kZ) = v{kZ, 0, 0) = - A — ' 1 4 ^ 

„ . .. w(kZ,0~b/2) + wae.O+fe/2) + w(kZ,-b/2 J)} + HjkZz+bj 2, 01 

w(kZ) = w(kZ, 0, 0) = 1 L 4 

.. 1 r w(fc£.+fc/2.0) - w(kZ,-bl2,0) _ v(fcg,0,+fe/2) - v(fc£ ,0,-fr/2)l 

9i(*£) = 2 L b b J 

■ u(kZ^+b!2) - n(fct,0,-fr/21 

92 '**' = b 

, ,, P v _ u(kZ,+b!2,0) - u(kZ-bt 2,0) 

Y3^J “ b 


where k is an integer and 0 < kZ < L. Allowing for m controllers at k = k\, k 2, k m , 
the corresponding relations can be represened by a 6mxl20 matrix acting on the state 
vector. (We consider in this paper m ■ 1, 2 or 3.) Let B denote the transpose of this 
matrix. Then the state space dynamics with co-located sensors can be described by: 


Mx + Ax + BU + BN a = 0 


(3.1) 


with sensor data: 

v p (f) - B*x(t) + N p (t) (3-2) 

\ r (t) = B*x(t) + N r (t) (3-3) 


where (/(•) denotes the control; N a ('), N p ( •), N r { •) model additive noise taken as 
(mutually independent) white Gaussian with spectral density matrices d a I, d p l , d r I 
respectively, / being the identity matrix. 
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4. The Continuum Model 

As we noted in Section 1, the problem of producing an “exact,” “three-D” continuum 
model for a real-world structure like the truss we are considering can be a formidable one 
although research in this area looks promising [5]. One way out of this difficulty is to 
exploit where possible the special nature of the truss — in our case it is a lattice of bays 
along the same axis numerous enough so that it is even visually “beam-like.” In that case 
there are many ways to approximate by “one-D” beams — without going into the details 
of this theory, suffice it to say that the approach by Noor and Russell [2] is the one 
adapted here. We thus create an “equivalent” (referring to [2] for the precise sense) one- 
dimensional anisotropic Timoshenko beam as follows. 

u denoting axial (longitudinal) displacement (along the x r axis); 

01 the torsion angle about this axis , 

w, 0 2 denoting the transverse bending displacement in the x r x 3 plane 
and torsion angle about the x 2 -axis , 

v, 0 3 denoting the transverse bending displacement in the x r x 2 plane 
and torsion angle about the x 3 -axis , 

the three axes being mutually perpendicular; 0 < x\ ^ L, 

L being the beam length = ni ; n = number of bays 

The Timoshenko equations (valid between control nodes) are: 

m u 'u - c n w" - c 14 v" - c 15W'” - Ci 502 + c 14 03 = 0 

m 2lV - C44V" - Ci 4 u" + C 44 03 = 0 
- C55W" - c 15 w" - C 55 02 = 0 

m 44 0i - c 66 0j' - C 36 02 - C 2 60 3 = 0 
m 5 5<t>2 + m seh + c ]S u' + c J5 w' - c 36 <t>;' + c S5 <|> 2 - c 33 < t>2 - C 23 <t> 3 
"*66<t>3 + m 56h ~ c 14«' - C44V' - C 26 <J>j' - C 23 <t> 2 + C 44 <t> 3 - C 22 <(> 3 

where the superdots denote time-derivatives and the primes, space derivatives. The coeffi- 


= 0 
= 0 
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cients of these dynamic equations are related to the truss parameters as follows: (cf. [2]). 
The mass coefficients are given by: 


mu 


m 44 


4m b + 4 m t + 4m^ + m 6 

m 2 2 = m 33 = £ 


= 2m55 * 2m^6 * 


i(Sm b + 12m e + 8 m d + m 5 ) 


6ji 2 


tm 6 

m 56 - 


The stiffness (flexibility) c t j are given by: 


where 


4 tS b S d \i 2 

Cn = + s d + s b (t + i?) 


c 44 - 


'14 


= C 55 = 


'15 


2 ts b s d 


s d + s b a + n 2 ) 


« 3 S t l 3 S b S d 

C22 - C33 = ^2 + 4{ s d +s b (t + \i 2 )) 


t-23 - 


4(S d + 5 t (£ + It 2 )) 


C66 = ”2C26 = -2C36 


* 3 S„S. 




H(S d + S b {Z + (i 2 )) 


t 

* - b 


In order not to complicate matters unduly in this demonstration, we shall freeze all 
parameters except the cross-sectional area A which will then be the structural parameter to 


be optimized. In this case 


(40 + 24^2 )EA 
C " 9 + 4^ 

(Newton) 

2 ea 

-C 15 - C 44 = C 55 - , + 2^2 

(Newton) 

(2725 + 1476V2)E4« 2 
C33 = 2628 + 1336VI 

(Newton)m 2 
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(97 + 140^2 )EAl 2 ' 
t2i ~ [2628 + 1336V2 , 

(Newton)m : 

-c !6 . ic . 33V2 )EAt 2 

2 296 + 130V2 

(Newton)m' 

m \\ = m 22 = = (8 + 5V2)Ap kg/m 


- 217,55 - . go ^ 9^Mt z P 

kg * m 

(AZ 2 p) 



Once the coefficients c i} , m i; are defined (on whatever basis), we can develop the 
generic state space dynamic model analogous to the FEM formulas (3.1), (3.2), (3.3): 

M'x + Ax + Bu + BN a =0 (4.1) 


yo = B*X([) + N p (t) 1 
v,(r) = B*x(t) + N a (t) j 


(4.2) 


where N a (-), N p (-), N r (-) are white Gaussian noise with spectral density d a I , d p I, d r I 
respectively. 1 Only, the dimension of the state x(t) is not finite. The technique of 
derivation is also different, in particular in the role of the energy. See [6, 7] for details. 
Here we can only summarize the basic results. 


Case 1: One Controller 

We begin with constructing the state space model for one controller (“midcontroller”), 
and offset masses at each end, referring to [4, 6] again for more details and to [7] where 


the general case of distributed control is treated. Thus the state x(t) at time t is defined 


by 


x(f) 


/(•>') 
m t) 
/(•* 2 > 0 
f(L, t) 


t See [7] for generalization to arbitrary diagonal matrices. 
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where s parametrizes the beam axis, zero denoting one end and L the other, L being the 
total beam length, and s 2 denotes the location of the mid-controllers and /(•, f) denotes a 
(6x1 vector) function of s, 0 < s < L, representing displacements and angles: 


/(S, 0 


u(s,t) 

v(s,r) 

w(5,r) 

<t> 2 (i,r) 

<t>3(-y,0 


The stiffness operator A is defined as follows: 


Ax = 


8 

v 


X - 


/(•) 

/( 0 ) 

/(•O 

m 


g(s) = -A 2 f"(s) + At/Xs) + Aof(s) , 0 < s < s 2 ; s 2 < s < L . 


The derivative /'(•) has a discontinuity at s = s 2 , and 


v 


-L i m-A 2 f(0) 
-A 2 (f\s 2 +) -f\s 2 -)) 
Lif(L) + A 2 f\L) 


and thus defined, the potential energy 


= [Ax, x) = 



/'(*) 
As) ’ 


where 


f'(s) 

As) 


ds > 0 


H 


Ci 0 

0 C 3 
0 0 
-C\ 0 


0 -C 2 
0 0 


A 2 


Cl 0 

0 c 3 
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C \\ 

C 14 

c 15 

Cx = 

C\A 

C44 

0 



0 

C 5 5 


j ^66 

^36 

C 26 

C 3 = 

c 36 C33 

c 23 


Cll 

c 23 

c 22 

/l, = 

0 c 2 

-Cf 0 


1 

0 -C15 

c \4 

II 

0 

0 

C44 

1 

0 -c 5J 

0 


B is defined by 


and 


^0 = Diag. f 0, 0, 0, 0, c 55 , c 44 ] 


L x 


0 -C 2 
0 0 


BU 



U is 6x1 


B*x = f(s 2 ) 

The mass/inertia operator M is defined 

I A*o/(0 I 


Mx 


M b ,om 

M c f(s 2 ) 

M biL f(L) 


(6xl) . 
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where 


M 0 *= 


m \\ 

* m 2 2 


™33 


m 44 

W55 

m56 ^66 


Diag. Af^ o = ( m 0 , m 0 , m 0 , Diag. /o ) 

where m 0 is the offset mass at j = 0 and 7 0 is its moment of inertia about zero, and 
similarly 

Diag -Mb L " ( m L> m L> Diag. / L ) . 

See also [6] for more on M bo , A/^. 

m c 0 0 0 0 0 

0 0 0 0 0 

0 0 m r 0 0 0 

M, = 

c 0 0 0 

0 0 0 / c 

0 0 0 


where m c is the force actuator (rotor) mass and l c the moment of inertia of the moment 
actuator about its center of gravity. The m.s. attitude error matrix is defined by 

1 im { i f /( 0, r)/(0, )* dt + \ J f(L, t)f(L, 0* dt\ 

T -*<*>{ * 0 0 ' 

and is the sum of the diagonal terms as defined. 


Case 2: Two Controllers 

Next we consider the case of two controllers, one at each end. Here, since there is 
no mid-controller, we may delete that entry in the state. Thus 


x(t) = 


/(■,0 
m 0 

f(L, t) 
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Mx 


M biL f(L,t) 


where the end-masses m 0 , m L must now include the actuator moving masses, and 
similarly for the moment of inertia matrices. We shall use the notation 


M bt0 0 

0 M b , L 


With U (- ) denoting the control vector, (12x1), we have 


Finally 


BU 


0 

U 


B*x 


/(0, f) 

f(L, t ) 



-L,/(0) -A 2 f(0) 
f(L) + A 2 f\L) 


8( s ) — —Aif + A\f’(s) + Aq/(s) , 0 < s < L . 

Here the mean square attitude error-matrix 

1 r T 

= lim (B*x(t))(B*x(t))* dt . 

T - »°° 1 n 


Case 3: Three Controllers 

In this case we have a mid-controller at s = s 2 as well as a controller at each end. 


Here 


x(t) 


/(•, t ) 

/( 0, t) 
f(S 2 , t ) ' 

f(L , f) 
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in other words the same state vector as in Case 1. 


Mx{t) = 


M 0 /(-, 0 
M b o /(0> 0 
M bt2 f(s 2 , t) 
M b>L f{L, t ) 


where 


M, 


0,2 


m c 0 


0 m c 
0 0 
0 0 
0 0 
0 0 


0 

0 

m c 

0 

0 

0 


0 

0 

0 


0 0 
0 0 
0 0 


where m c is the force actuator moving mass and I c is the moment of inertia about its 
center of gravity, corresponding to the “mid-controller.” 


BU = 


B*x = 


0 

U 


AO, t) 
f(* 2- t) 
f(L , t ) 


We can calculate from the diagonal terms of the mean square error matrix: 

. T 

Diag. lim ^ J (B*x(t))(B*x(t))* dt . 

T - *» 1 o 

which now has 18 entries. Once again we adopt the notation: 


M n 


M r = 


Mi 


h,0 


M 


42 


M 


b,L 
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5. Formulas for Performance Criteria 

We shall now develop formulas for the Performance Critera. First the mean square 
attitude error: Using either model FEM (3.1) >(3.3) or Continuum (4.1)-(4.2) we have the 
state space form: 

M'x + Ax + BU + BN a = 0 

v p - B*x + N p 

\ r * B*x + N r 

and 

- CLWp(t) + yv r (0 . 

Substituting this control law into the state equations we have: 

Mx + (/! + a£fl*).x + y££*j; + B(A/ a + aA^ + yN r ) « 0 . (5.5) 

The steady-state output covariance matrix 

[ o?d + d +y^d a 

2 2^ + aBB*)~ l B) (5.6) 

where the matrix part: 

B*(A + <xBB*)~ l B 

is recognized as the steady-state input-output response matrix: it is the value at co = 0 of 
the input-output transfer-function: 


(5.1) 

(5.2) 

(5.3) 

(5.4) 


B*((id)) 2 M + A + a BB* + (i(ti)yBB*)~ l B . 

The scalar factor 

J + d a +y 2d r 

d ' l-i < 5 . 

consolidates the “noise” part. For the case of two controllers the mean square attitude 
error can be calculated from 

Diag. R a . 

For the third case of three controllers it is given by 
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Diag. LR a L* 


where 


LU 


m 

m 


For Case 1 with one controller only, we have to calculate 


£(/( 0 , 0 /( 0 , 0 *) 


£(/(£. 0 *) • 


Separately, expressing each as a transformation of the state. 

L\x - /( 0 ) 

Lix - f(L) . 

In the FEM version, we have thus to invert the matrix 

(, A + a BB*) 

which in our case is 120x120 — and of course can be done only numerically. For the 
continuum model however we have to invent the operator 

(A + clBB*)~' 

but — and this is the main point of departure — this can be done analytically. Referring 
to [7] for details, here we shall simply enumerate the formulas below. 


Case 1: One Controller Only 

_. f 1 l+.?i ln-rl 1 1 i]j 

E[/(0, r)/(0, r)*] = Diag. a ’ a ’ a’ a’ a) a 

,, fl 1 + (L-si) 2 1 + (L-Si) 1 1 1 L) d 

E[f(L,t)f(L,t )+ ] = Diag. a - a ’ a ’ a ’ aj 

(each is a diagonal matrix!). Steady state step response corresponding to step input U: 
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m~) = £ 


0 

0 

0 


0 0 
0 -j 2 

S 2 0 


u 


AL, «) - i 


Note: Controller at j = j 2 . 


0 

0 

0 


0 0 
0 L-j 2 
s 2 -L 0 


U 


Case 2: Two End Controllers 

Diag. £[/(0, r)/(0, r)*] = Diag. £[/(L, t)f(L, r)*] . 
And the first four diagonal terms in either matrix are given in order by: 

-[>♦$$] 


where 


c 33 


+ vL . £*&_ 

t /«| “T A 


2a + 4 


( c » *f + ! <‘)( c “ + ¥)- c -? 


4 = 


aL 


2c6 4 

P - 4 + 1 

X = n + 2 + ( g n ~ 2c)(2a + cSL) 

2ac8 
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c 


— C44 


6 = 


7 = 


1 + 


Hr 

12 


£m_ 

C33 - 


^66 


Case 3: Three Controllers 

Here we have to express the answers in terms of 6x6 matrices. Thus let 


D „ 

D 2 2 


D 


33 . 


£[/(o, 0 m 0*] 

d 

£[/(*2, 0 /(*2, 0*] 

d 

£[/(L, r) /(L, <)*] 
d 


Then 


Du 
D 22 
D 33 


= ^llfo) _ <fl2(^2)(^22(^2) + (dll(L-S 2 ) 1 “ tt) ) ^2(^2) 


(^22(^2) + d\ y{L-s 2 ) 1 - cx) 

(a/ + m 2 {L-s 2 ))~' {a/ + m 3 (L-s 2 ) + m 2 l (L-s 2 )(D 22 )mt\(L-s 2 )} 

X (a/ + m3(L-J2)) 1 

<*ii(s) - ((a + m^s)) - m 2 l (s)*(a + miis ))' 1 m 2 i(s))~ l 
d 22 (s) = ((a + m 3 (s)) - m 21 (s)(a + m 1 (j)) _1 mJi(.r))~ 1 


N-l 


where 


d n (s) = -d u (s)m 2] (s)*(al + ntiis)) 


-1 


mjCr) = 




c5fc) 


M2 


cSQ) A * 

2 


An 


C 3 cj 5(5) 4 

+ — 7^ A 2 2 
s 4 


m 3 (s) = 


5 

£§tel A 
2 


1 1 

- 

C % S ) A t - 
2 An 

* 

12 

S 

cs5(s) . 
~ 4 A 22 
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where 


- f A*,, 

m 2I (s) = 

c5(r) ,* Cj cs8(s) A 

I 2 A ' 2 "^ + 4 ^ 



0 0 0 

A 22 =010 
0 0 1 


Stability Margin Formula 

From (5.5) it is clear that the closed-loop eigenfunctions are specified by 
\ 2 M$ + (A + a£B*)§ + y\BB*$ = 0 

(again, irrespective of whether we are using the FEM version (3.1) -(3.3) or the 
continuum version (4.1), (4.2)) where <)> denotes the eigenvector (“mode shape”), with the 
eigenvalue X specified by 


X 2 [A/<t>, <|>] + [(A + a BB*)Q>, <|>] + = 0 

where we may normalize the mode shape <|> so that [M<J>, <M = 1 . Since a 2: 0 and y > 0, 
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it follows that the eigenvalue X is given by: 


where 


X - -a ± ico 

a _ v Jg!M i 

o - -T 2 
Cl ) 2 = COq - a 2 

= [(^ + . 


Thus the closed-loop eigenvalues have strictly negative real parts. 

It is shown in [8] that the sum of the absolute values of the closed-loop eigenvalues is 

given by 

Y Tr. M; 1 (5*8) 


in all cases, where in M c only moving parts of the actuators must be considered (as 
opposed to the stationary mass such as the armature mass). The simplicity of this formula 
is striking when compared with taking the sum of the inverses of the absolute values of 
the real parts of the closed loop eigenvalues for a finite number of modes as in [10]. We 
may note that for zero natural damping (or with damping, if we consider only the increment), 



for any continuum model. Again, (5.8) applies for the FEM version (3.1) -(3.3), as well 
as for the continuum model (4.1), (4.2). 
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6. Numerical Results 


In this section we compare the numerical results by the two methods for the mean 
square attitude errors (and equivalently, tip response), for each of the three cases: one, 

two and three controllers — using the nominal values (cf. Section 1) for the truss para- 
meters. 


Case 1: One Controller 
FEM 

The tip response /( 0, °°) and f(L, °°) was calculated for: 


c 44 = 35571851.2 ; 

& 

II 

O 

L L L 
9 ’ 3 ’ 2 ’ 

10,000 ; 

Ni 

II 

O 

L L L 
9 ’ 3 ’ 2 ’ 

100 ; 

O 

II 

<N 

L L L 

9 ’ 3 ’ 2 ’ 


both /( 0, -) and f(L, «) values were exactly the same as the values predicted by the 
continuum model and hence are not displayed. 


Case 2: Two Controllers 


*« - Diag. »'» + alf' B , , . i 12 




FEM 

Continuum 

a = 1 

b x 1 = b-n 
&22 = ^33 

= 6gg = b 9 g 

b* 4 = ^ 10,10 

0.4998 

0.99727 

0.50028 

0.5 

0.99727 

0.5 

a = 50 

b\ ] = b 11 
&22 * £33 

“ ^88 = ^99 
b 44 = b]QiQ 

0.01 

0.019945 

0.01 

0.01 

0.019945 

0.01 

CC = C44 

= 3.557xl0 7 

b } 1 = b-n 

&22 - ^33 

~ ^88 “ bgg 

b 44 = ^ioio 

2.2978 xl0~* 
2.80545x10'* 
2.478 x 10 * 

2.2978 xlO'* 
2.80546x10'* 
2.478 xlO'* 


27 Q 



Case 3: Three Controllers 


a - 3.557 x 10 7 
b u = « = 1 , .... 18 


Actuator 

Position 

b n x 10 8 

oo 

o 

X 

FEM 

Continuum 

L 

s ~ 9 

1.5119 

1.5127 

3L 

Sm T 

1.8366 

1.8370 

AL 

s " T 

1.9461 

1.9465 

II 

2.0330 

2.0333 

II 

2.1031 

2.1033 

1L 

S ~ 9 

2.1602 

2.1603 


In other words the FEM and the continuum gave exactly the same numerical results 
within (the SUN-386i) computer accuracy in all cases. 
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FIGURE 3 

MIDACTUATOR POSITION: INCREMENTS L/16 
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SUM OF M.S. BENDING DISPLACEMENTS AT BOTH ENDS 



FIGURE 4 

MIDACTUATOR POSITION: INCREMENTS L/16 
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SUM OF M.S. TORSION AT BOTH ENDS 



FIGURE 5 

MIDACTUATOR POSITION: INCREMENTS L/16 
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M.S. AXIAL DISPLACEMENT AT ACTUATOR LOCATION 



FIGURE 6 

MIDACTUATOR POSITION: INCREMENTS L/16 
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7. Optimal Location of Controllers 

As a byproduct of our theory, we can examine the problem of the optimal actuators 
placement that minimizes the mean square attitude for any given choice of control gains 
and structure parameters. 

Case 1: One Actuator 

Here we take the criterion as the sum of the mean square displacements at both ends. 
We calculate explicitly that 

[£[/(0,f) /(0,r)*] + E[f(L,t) f(L,t)*]] 

d 

\2 2 + s 2 + (L-g) 2 2 + s^ + (L-sj 1 2 2 2] 

- Diag. |_ a , a . a ’ a ’ a ’ aJ 

from which it is clear that the optimal placement is in the middle: 

L 

s ~ 2 ’ 

Only the bending-displacement is affected by actuator position. 

For Case 2 of two actuators there is no placement problem since one is required to 
be at each end. 

Case 3: Three Actuators 

Here we can consider the optimal placement of the midactuator while the other two 
are fixed one at each end. The behavior of the sum of the mean square errors (setting 
a = 1) at both ends for the axial mode is shown in Figure 3; and the bending mode is 
shown in Figure 4; and the axial torsion in Figure 5. In all cases we see that the worst 
position is at the middle! The best place is at either end. Finally Figure 6 shows the mean 
square axial displacement at the actuator location. Again the worst place is the middle. 
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8. Optimization 

We shall now treat a canonical optimization problem currently studied by FEM and 
truncated modal models [9, 10]. The objective is to minimize the total mass — structure 
and control — subject to meeting specified indices of performance. Here we take them to 
be: 

(a) mean square attitude error due to sensor noise less than or equal to fixed value 

(b) stability index: sum of the absolute values of the real parts of the closed loop 
eigenvalues to be not less than a fixed value. 

We shall see that the problem can be solved analytically by virtue of the formulas we 
have developed using continuum models. 

The structural parameter we shall use is the cross-sectional area A of the longerons 
(assumed to be the same for battens and cross-bars). Other parameters being fixed, the 
structure mass is then proportional to A. (The extension to the case of nonequal areas 
only complicates the algebra, as can be seen from the expressions (cf. Section 3) for the 
flexibility coefficients.) The control mass has to be subdivided into a stationary mass 
(armature mass, for example) and a moving mass (rotor mass, for example) since only the 
latter is involved in the stability index formula. The stationary mass is of course related to 
the moving mass — for simplicity we shall take it to be inversely proportional to the 
rotor mass. The control parameters are the attitude and rate gains a and y. These of 
course will need to be constrained not to exceed prescribed limits. Thus we have the 
following formulation (nominal values for all structure parameters except A): 

Structure mass = NZpA 

hr 

Control stationary mass = ^ 

Moving mass = m 

h 

Total mass *= NZpA + — + m 
r m 

Stability index = — 

J M + m 
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(where M denotes the contribution of the end-masses). For the truss considered, 


N = (76 + 46^2) . 


Finally the mean square attitude error — to be specific, we shall consider the case of 
two controllers, one at each end; and take the sum of the mean square displacements at 
either end. First we express these in terms of the structural parameter A — we have thus 
to use the expressions we have derived for the flexibility coefficients {c tJ } in terms of A 
in Section 5 and substitute them into the formulas for mean square errors for two controllers 
in Table 3 under Case 2. In doing so we shall also take advantage of the simplification 
possible by noting that for the nominal value of L = 27 meters, we can readily calculate 
that 

Xp > 2 


so that we may replace both (X.p — 1) and (A.p - 2) by X.p. Thus the first four diagonal 
terms in 

2 Diag. £[/( 0, r)/(0, 0*1 
( = 2 Diag. £[/(L, f) /(£, 0*] ) 

are given in order by: 


/(<*, y, A) = 


a 


1 + 


U_ L a 
4 + 2c5 EA 




a 


X 

a 


1 + 


L a_ 
2 c8 EA 


1 + 


L _a_ 
2c5 EA 


X 

a 


1 + 


L _a_ 
2c8 EA 


1 + 


L _CL 
2cS EA 


a + 


c 33 + 


c 66 


L 

2 EA 


2 EA f c 66 L au 

[C33 + -f + 2 J 1 


L a 3 


C66 + 2 £/lJ 


256 


— ' 
2 J 



where now 


c \ i = 


40 + 24^2 
9 + 4V2 


2 

1 + 2V2 


2725 + 1476V2 
2628 + 1336V2 


32 + 66^2 
296 + 130V2 


n = number of bays ; L = nZ 


12 (c 33 - c 66 / 4) 


with <i as given by (5.7). These formulas enable us to draw conclusions concerning the 
dependence on the cross-sectional area A without resorting to numerical computer calcu- 
lations. We see that all the errors decrease as A increases. The axial error decreases from 
2d d 

— at A = 0 to — at A = <*>; similarly the torsion error. The bending error is least 
affected, decreasing from ^ at A = 0 to j_7 j ~ at A = in all cases the 

minimal mean square error is at most 3 db less than the maximum! 

For the optimization let us fix on the mean square bending error as being the easiest 


analytically: let 


d 1 

/(a, y. A) - - 1 - -n = 

a Lf_ L a 

4 + 2c8 EA. 


Thus the optimization problem is that of minimizing 


subject to: 


NZp A + ^ + m 


M + m 


/(a, y. A) £ cl 
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The first inequality can clearly be reversed to read 


M + m 
Y 


< 


1 




We note that the objective functional 

NlpA + £ + m 

is infinitely smooth and trivially convex, and the constraints are also infinitely smooth 
and convex. Hence we are assured of the existence of a minimum (which is further 
verified to be unique). Moreover we can go to the Langrange parameter formulation and 
minimize: 


(/V£p.4 + ^ + m] + Xj J + ^/(a, y, A) 

where t X 2 ^ 0 are the Lagrange parameters. See [11] for the standard results that 
are applicable here. 

Compared to the FEM versions [9, 10] this is a “trivial” problem and complete 
“analytical” solution is possible. We omit the details since our primary aim in this paper 
is to demonstrate the simplicity of the optimization problem relative to the FEM versions. 
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